Calculated Magnetic and Electronic Properties of Pyrochlore Iridates 
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Using density functional theory and LDA+U method, we investigate magnetic and electronic 
structure of Y2lr207 and rare-earth based pyrochlore iridates. Our study reveals that the ground 
state is a non-coUinear magnetic insulating state. Due to strong spin-orbit coupling in Ir 5d, there is 
an unusual correlation between the bands near Fermi level and the magnetization direction, resulting 
in a possibility of insulator-to-metal transition under applied magnetic field. This makes pyrochlore 
iridates a good candidate for possible magnetoressitance and magnetooptical applications. 
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While it has been realized that the Coulomb in- 
teraction is of substantial importance in Srf transi- 
tion metal systems, resulting in a great variety of 
their physical properties, such, e.g., as metal-insulator 
transition^], colossal magnetoresistance[2j and high-Tc 
superconductivity 0, the 4d and the 5c! orbitals are, on 
the other hand, spatially more extended and usually pos- 
sess a much broader bandwidth, making the correlation 
effects to be minimal. However, very recently, both the- 
ory and experiment give the evidence on the importance 
of interplay between spin-orbit coupling and Coulomb 
interaction for several, mainly Ir based, bd compounds, 
which have received a lot of research attention |4l4l4 1 . In 
particular, in Sr2lr04, one of most studied Ir oxides, it 
has been found that the material becomes a Mott insu- 
lator 0, @|, and exhibits very interesting magnetoelec- 
tric properties^. A different system, Na4lr308, which 
crystallizes in a geometrically frustrated spinel structure, 
is an insulator with a large Curie-Wiess temperature 
(about 650 K) and with a considerable effective moment 
(1.96 ^B)[i5[. However, this compound does not exhibit 
any sign of rnagnetic order even for the lowest measured 
temperature [l5|, and it has been suggested as one of the 
few long sought quantum spin liquids [B-lli. There is 
also theoretical work that addresses anisotropy of mag- 
netic interactions [igl. [20| . 



out structural change [2l|. There are also several elec- 
tronic structure calculations trying to explore electronic 
and magnetic properties of the iridates with geometri- 



Ir oxides ^2lr207|6!, |21|-|26|| {A^ Y or rare-earth el- 
ement), which crystallize in a pyrochlore structure 27| , 
is another geometrically frustrated iridate system. It 
has been recently discussed in connection with a novel 
"topological Mott insulator" phase [ij seen between 
topological band insulatingfisi] and conventional Mott 
insulatingjj phases as interactions get stronger. Exper- 
iment observes that depending on the A^iie . A o h-207 
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show a wide range of electrical properties 
example, Y2lr207 is an insulator^fll but with increasing 
the ionic radius at the ^-site, the system eventually be- 
comes metallic for Nd2lr9 07 2l|. while Pr2lr207 shows 
strong Kondo behavior |25l|. Moreover, it has been found 
that temperature will drive an insulator-to-metal tran- 
sition associated with abnormal magnetic behavior with- 



cally frustrated structure [19|, 12 



In the present work, we perform a detailed study of 
the magnetic and electronic structure for ^2li'207. We 
find that the ground state of those systems is a non- 
coUinear magnetic state while geometrically frustrated 
pyrochlore lattice makes other states with different ori- 
entations of moments close in energy. When ^=Pr and 
Nd the electronic structure of these compounds shows 
metallic behavior while for A—Y, Sm, Eu the materials 
are insulators. We uncover that spin-orbit (SO) coupling 
affects the energy bands near the Fermi level and depend- 
ing on orientation of moments some of those systems can 
be switched from an insulator to a metal. Exotic elec- 
tronic and magnetic properties would make Ir pyrochlore 
to be good candidates for various applications including 
magnetoresistance effect and magnetooptics. 

We perform our electronic structure calculations based 
on local spin density approximation (LSDA) to density 
functional theory (DFT) with the full-potential, all- 
electron, linear-muffin-tin-orbital (LMTO) method (sot. 
We use LSDA-I-U scheme [sij to take into account the 
electron-electron interaction between Ir 5d electrons, 
and use U = 2 eV which has been previously found to 
be adequate in iridates [3, M, 0- When the A site is a 
rare earth element, we also add the Coulomb interac- 
tion for the localized 4/ electrons and use U = 6 eV. 
We use a 12x12x12 k-mesh to perform Brillouin zone 
integration, and switch off symmetry operations in order 
to minimize possible numerical errors in studies of vari- 
ous (non-)collinear configurations. We use experimental 
lattice parameters [i^l in all set ups. 

We first discuss our results for Y2lr207. Without spin 
orbit (SO) couphng, both LSDA and LSDA-I-U predict 
this system to be metallic which is not consistent with 
the experiment. Since the SO strength is large for Ir 5rf 
electrons (about 0.4 eV)[32|, and has been found to pro- 
duce an insulating behavior in Sr2lr04[3,S|, we have per- 
formed the LSDA-f-U-t-SO calculations. There are four Ir 
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atoms inside the unit cell forming a tetrahedral network 
as shown in Fig.l. We first set an initial magnetization 
axis along (001) direction, however the calculations con- 
verge to a non-collinear state with the magnetic moment 
departing from the initial orientation a little bit. The 
four Ir sites have similar spin and orbital moments (about 
0.13 fig) and both of them are slightly smaller than the- 
oretical values reported previously^]. However, due to 
the obtained slight non-collinearity of the solution, the 
net magnetic moment is found to be small. 

Each of four Ir atoms is octahedrally coordinated by 
six O atoms, which makes the Ir 5d state split into doubly 
degenerate e^ and triply degenerate t2g states. Due to the 
extended nature of Ir 5d orbital, the crystal-field split- 
ting between t2g and Cg is large with the e,, band to be 2 
eV higher than the Fermi level. The bands near the Fermi 
level are mainly contributed by Ir t2g with some mixing 
with O 2p states. SO coupling has a considerable effect 
on these t2g states: it lifts their degeneracy and produces 
24 separate bands in the range from -2.3 to 0.7 eV. Same 
as in Sr2lr04[l, the bandwidth of these t2g states in 
our LSDA+U+SO calculation for Y2lr207 is also narrow, 
however, it is still metallic as shown in Fig. 2(a). Naively 
one may expect that using larger Coulomb U will result 
in an insulating state. However, our additional calcula- 
tions show that increasing U cannot solve this problem, 
and even a quite large U (=5 eV) cannot open a band 
gap for the initial coUinear (001) setup. Here we agree 
with the previous calculation [29j, which also showed that 
U cannot open a band gap. 

Strong spin-orbit coupling may however induce the 
dependence of the band structure on quantization axis. 
So, we subsequently perform our calculations with initial 
magnetization aligning along (110), (120), (111) direc- 
tions. We also perform the calculations with two sites in 
a tetrahedron along and other two pointed oppositely to 
(001), (111), (110) or (120) direction in order to account 
for possible antiferromagnetism. Interestingly, we find 
that despite the strong SO coupling which usually re- 
sults in a large magnetic anisotropy energy, for Y2lr207, 
the rotation of magnetization does not involve too much 
change in the total energy. Among the magnetic configu- 
rations mentioned above, the (111) direction is found to 
be lowest, but the energy difference between this and the 
highest energy (001) state is just about 3.7 meV per unit 
cell. It is also interesting that while the magnetization 
direction does have a considerable effect on the bands 
near the Fermi level, as can be seen from comparing the 
Fig. 2(a) and Fig. 2(b), it turns out that all of the above 
mentioned calculations converge to metallic states. Also, 
all of them produce a considerable net magnetic moment 
in contrast to the experiment 22|424l|. Thus, one can con- 
clude that neither of these magnetic configurations may 
be the ground state. 

With the pyrochlore structure, the Ir sublattice has a 
topology consisting of corner-sharing tetraliedra and is 



TABLE I: The spin {S) and orbital (O) moment (in fig), the 
total energy Etot (in meV) as well as the band gap Egap(in 
meV) for several selected magnetic configurations of Y2lr207. 



Configuration: 


(001) (111) 2- 


-in/2-out all 


-in/out 


(S) 


0.13 0.17 


0.14 


0.13 


(0) 


0.13 0.13 


0.15 


0.17 


Etot(meV) 


4.47 0.69 


3.21 


0.00 


Egap(meV) 
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geometrically frustrated. Thus, we carry out several non- 
collinear calculations with the initial state to be "all- 
in/out" (where all moments point to or away from the 
centers of the tetrahedron), "2-in/2-out" (two moments 
in a tetrahedron point to the center of this tetrahedron, 
while the other two moments point away from the cen- 
ter, i.e. the spin-ice [33| configuration), and "3-in/l-out" 
magnetic structures (see Fig.l for the moments configu- 
ration). We find that the "all-in/out" configuration is 
the ground state, and this state is insulating as shown 
in Fig. 2(c): it has a band gap of 0.03 eV. The exper- 
imental band gap for Y2lr2 07 is not available, though 
it had been expected that it should be larger than that 
of Sm2lr207 (0.01 eV)'2l']. Different from other mag- 
netic configurations, during the self-consistency the "all- 
in/out" state will retain their initial input direction; thus, 
there is no net magnetic moment. This is consistent with 
the experimental fact on the absence of the magnetic 
hysteresis loopjlJ]. As shown in Table I, the energy dif- 
ference between the ground and several selected excited 
states with different orientations of moments is small. 
Here, the energy difference involves not only the energy 
of magnetic anisotropy but also the effect of intersite ex- 
change interaction. Due to the geometrically frustrated 
structure this energy difference is obviously smaller than 
that found for another 5c? compound Ba2NaOs06[l2!|. 

Regardless the proximity of the ground and excited (or 
said metastable) states in energy, they would have very 
different conductivity and magnetic properties. Thus one 
can understand that for the same compound Y2lr207 
Taira et aZ.|24| observe no ferromagnetic ordering while 
Yanagishima et al. 22 , [i^ claim a presence of a small 
net magnetic moment. One can also understand the ob- 
served temperature induced metal-insulator transition, 
as well as a large difference in temperature dependence 
of magnetization measured under zero-field-cooled con- 
ditions (ZFC) and under field-cooled conditions (FC) at 
low T's. 

Based on both strong sensitivity of the energy bands 
near the Fermi level on the orientation of moments, and 
proximity of various magnetic states in energy, it is natu- 
ral to expect that an application of a magnetic field could 
have a big effect not only on the magnetic response but 
also on the conductivity in iridates. In particular, this 
should result in a large magnetoresistance effect if one is 



FIG. 1: The pyrochlore structure showing Ir tetrahedral net- 
work and its magnetic configurations (a) "aU-in/out" config- 
uration, (b) "2-in/2-out" configuration. The arrows denote 
moments directions. 



able to switch between insulating "all-in/out" state and 
any collinear state. This simple idea has been proved 
by the following numerical calculation. Starting from 
the "all-in/out" ground state, we apply an external field 
along (001) direction. The result shows that the exter- 
nal field will rotate the magnetic moments meanwhile 
only slightly change their magnitude. AST magnetic 
field along (001) induces a 0.07 /Xg net magnetic mo- 
ment, which is in fact close to the experiment performed 
for Sm2lr207, where it was shown that a 4 T magnetic 
field produces a 0.05 total moment [2^ . Increasing 
the field further, the numerical calculation does find en- 
ergy bands crossing the Fermi level, namely an insulator- 
to-metal transition at a field of 40 T, although there is 
already non-negligible density of states at E/ for a lower 
magnetic field. 

Another interesting feature which emerges from this 
electronic structure is an expected magnetic field depen- 
dent optical response. We have computed interband op- 
tical conductivity crmter(w) of Y2lr207 for (001), (111) 
and "all-in/out" orientations of moments. The results 
are shown on Fig. 3 where one can see a very different be- 
havior of this function depending on the imposed mag- 
netic configuration at infrared frequencies. The intra- 
band Drude contribution can also be restored from our 
computed intraband plasma frequencies for the metallic 
(001) and (111) states which are pretty small and equal 
to 0.7 and 0.2 eV, respectively. Thus, controlling the 
strength and direction of the applied magnetic field one 
can reach a continuous change from metallic to insulating 
response in this system. 

Experiment finds that changing the 4-site of ^2lr207 
will vary the properties considerably. Yanagishima and 
Maeno 22[ show that for ^2li'2 07, the electrical conduc- 
tivities of ^=Pr, Nd, Sm and Eu exhibit metallic be- 
havior, while the electrical conductivities of ^=Gd, Tb, 
Dy, Ho, Yb and Y are like in insulators. Matsuhira et 
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FIG. 2: LSDA-I-U-I-SO electronic band structure of Y2lr207 
for difi^erent orientations of magnetic moments: (a) along 
(001) direction; (b) along (111) magnetization direction; (c) 
"all-in/out" configuration. 
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Ho, Yb and Y are like in insulators. 
a^.[2lj also observe that the A-s\te induces the change 
from metal to insulator, but they claim that it is taken 
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FIG. 3: Calculated interband optical conductivity of Y2lr207 
for different orientations of magnetic moments: black - along 
(001) direction; red - along (111) magnetization direction; 
blue - "all-in/out" configuration. 
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place around Nd. Namely, the systems with ^=Nd, Sm 
and Eu are insulators while those with A=Pr are metals. 



Matsuhira et aL[2l| try to contribute this discrepancy to 
the quality of samples. However, even using the same 
synthesis condition, the discrepancy still appears pi). 

To clarify the effect of ^-site, we first perform a con- 
strained calculation with the 4/ band shifted by a con- 
strained potential. This however almost does not affect 
the bands around the Fermi level. So, one can conclude 
that the rare earth element has a small effect on the 
conductivity. We further perform the calculation by us- 
ing ^2lr207 structure but replacing the rare earth ele- 
ment A by Y. Like in Y2lr2 07, the ground state is found 
to be "all-in/out" non-collinear solution. However for 
Pr2lr207 and Nd2lr207, the ground state is metallic. 
Changing the ^-site from Nd to Sm and further to Eu re- 
sults in decrease in the ionic radius. This reduces the Ir~ 
0-Ir angle, and makes our calculation for ^=Sm and Eu 
to produce insulating band structures. We contribute the 
discrepancy between Ref. 22| and Ref. 2l| to the small- 
ness of the energy difference between the "all-in/out" 
ground state and other metastable states. 

In summary, using the LSDA-f-U-f-SO method, we have 
explored the electronic and magnetic properties of ge- 
ometrically frustrated pyrochlore iridates. Our study 
reveals that the ground state of these systems has a 
non-collinear "all-in/out" magnetic configuration of mo- 
ments. Our constrained calculation shows that the A site 
has a small affect on the energy bands near the Fermi 
level while different ionic radii and corresponding changes 
in lattice parameters are the reason why ^=Pr are metal- 
lic while ^=Nd, Sm, Eu, and Y are insulating. Thanks to 
the strong SO coupling in Ir 5 c? states, rotation of magne- 
tization can produce insulator-to-metal transition thus 
making iridates to be interesting candidates for various 
applications, such, e.g., as giant magnetoresistance effect 
and magnetooptics. 
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